bST0-'m-0435’ 


Broadband  Active  Sonar: 
Implications  and  Constraints 

H.  Lev/ 


19970307  106 


Broadband  Active  Sonar: 
Implications  and  Constraints 


H.  Lew 

Maritime  Operations  Division 
Aeronautical  and  Maritime  Research  Laboratory 

DSTO-TR-0435 


ABSTRACT 


The  broadband  sonar  concept  has  an  impact  on  all  areas  of  an  active  sonar  system. 
This  includes  the  signal  processing  of  transmit-receive  signals,  the  response  of  the 
medium  and  the  target  response.  Beginning  with  the  definition  of  a  broadband  signal, 
the  implications  of  broadband  signal  processing  in  terms  of  matched  filtering,  pulse 
compression,  low  probabiKty  of  interception  and  beamforming  are  examined.  This  is 
followed  by  an  investigation  of  the  constraints  imposed  on  the  propagation  of 
broadband  signals  by  the  medium.  Finally,  a  brief  discussion  on  the  use  of  broadband 
techniques  for  target  classification  concludes  the  report. 
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Executive  Summary 

The  basic  motivation  for  using  broadband  sonar  is  to  improve  the  range  resolution 
capabilities  without  sacrificing  the  detection  performance,  which  is  governed  by  die 
amotmt  of  received  energy  from  the  target.  This  can  be  achieved  by  mcreasing  the 
bandwidth  of  the  signals  used.  From  this  improvement  in  range  resolution,  a  number 
of  other  desirable  features  follow  such  as  improved  range  accuracy,  robust  parameter 
estimation,  Doppler  tolerance  and  so  on. 

The  report  begins  by  examining  the  definition  of  a  broadband  signal.  The  fractional 
bandwidth  and  time-bandwidth  product  definitions  are  discussed.  This  naturally 
leads  on  to  a  discussion  of  broadband  signal  processing  which  includes  matched 
filtering,  pulse  compression  techniques,  low  probability  of  interception  signals  and 
beamforming. 

Having  examined  the  signal  processing  aspects,  it  is  then  necessary  to  determine  the 
response  of  the  medium  to  the  transmission  of  broadband  signals.  By  using 
propagation  models  or  experimental  data  the  usable  broadband  frequency  range  can 
be  determined  for  a  given  sonar  system  -  environment  scenario.  For  one  example  it 
was  foimd  that  a  spread  of  frequencies  between  100  to  2000  Hz  could  be  transmitted 
without  significant  distortion  across  the  spectrum  over  a  minimum  range  of  10  km 
given  a  probability  of  detection  of  90%  and  false  alarm  probability  of  one  in  ten 
^ousand. 

Finally  the  report  gives  a  brief  discussion  on  tiae  implications  of  target  classification 
using  broadband  techniques  such  as  resonance  scattering  and  high  resolution  range 
profiling. 
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1.  Introduction 


In  active  sonar  systems  pulses  with  a  single  carrier  frequency  are  often  employed 
because  of  their  simplicity  and  ease  of  implementation.  They  are  sometimes  referred  to 
as  "simple"  signals  and  have  the  property  that  their  time-bandwidth  product  is  of  the 
order  of  unity.  This  property  makes  long  range  detection  incompatible  with  high  range 
resolution  as  detection  essentially  depends  on  the  amount  of  energy  received  from  the 
desired  target.  Improved  detection  can  be  expected  by  increasing  the  peak  power 
and/or  the  pulse  length.  In  particular,  long  range  detection  requires  the  target  echo  to 
have  a  sxifficiently  large  amplitude  since  the  amplitude  decreases  with  range[l].  In 
practice,  however,  transmitter  limitations  (e.g.,  transducer  dynamics,  cavitation, 
interaction  effects  etc)  prevent  excessive  peak  power  outputs.  Alternatively  the  amotmt 
of  detectable  energy  can  be  increased  by  lengthening  the  pulse  duration  but  this 
comes  at  the  expense  of  reduced  range  resolution.  This  is  because  the  range  resolution 
is  determined  by  the  bandwidth  of  the  signal,  which  for  "simple"  signals  goes  as  the 
reciprocal  of  the  pulse  length.  To  avoid  this  inherent  conflict  between  the  needs  of 
reliable  detection  and  good  range  resolution  signals  other  than  the  "simple"  ones  have 
to  be  used.  These  are  signals  with  large  time-bandwidth  products  and  are  generally 
classed  as  broadband  signals.  In  this  report  an  overview  of  broadband  signals  and 
sy steins  in  the  context  of  active  sonar  will  be  given. 

There  are  several  areas  of  an  active  sonar  system  in  which  the  use  of  broadband  signals 
win  have  a  significant  influence.  For  the  purposes  of  analysis  and  discussion  it  is 
useful  to  divide  the  sonar  system  into  three  broad  subsystems  as  shown  in  Fig.  1. 


Figure  1.  Generic  active  sonar  system.  Note  that  the  transmitter  and  receiver  need  not  be  in 
the  same  location. 


The  transmitter-receiver  block  will  be  discussed  in  sections  2  and  3.  This  will  deal  with 
the  definition  of  broadband  signals  and  the  signal  processing  aspects  of  broadband. 
Section  4  will  look  at  the  impact  of  the  medimn  on  the  propagation  of  broadband 
signals.  In  section  5  the  implications  of  broadband  for  the  target  will  be  discussed  in 
terms  of  target  strength  measurements  and  classification.  The  relationship  between  all 
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three  of  these  areas  from  a  systems  viewpoint  is  tied  together  in  a  conceptual 
framework  shown  in  Fig.  2. 


Figure  2.  Conceptual  framework  of  sonar  system  modelling. 


2.  Broadband  Signals 

2.1  Fractional  bandwidth  definition 

Intuitively,  a  broadband  signal  is  one  which  has  a  significant  amoxmt  of  its  energy 
distributed  over  a  wide  range  of  frequencies.  At  present  there  is  no  standardised 
definition  of  what  a  broadband  signal  is.  The  terms  broadband  and  narrowband  can 
have  several  meanings  depending  on  the  application.  Note  that  broadband  is 
sometimes  referred  to  as  wideband  in  the  literature.  Other  terms  such  as  imptdse, 
nonsinusoidal,  baseband,  video  pulse  and  time  domain  are  also  used  in  connection 
with  broadband  signals  depending  on  the  application.  In  this  report  two  definitions  of 
what  a  broadband  signal  is  will  be  adopted.  The  first  is  a  quantification  of  the  intuitive 
notion  of  a  broadband  signal  and  the  second  is  derived  from  the  limitations  of  the 
sonar  signal  model. 
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Firstly,  consider  the  definition  of  the  bandwidth  of  a  signal.  The  bandwidth  is  the 
frequency  range  within  which  a  significant  fraction  of  the  total  energy  of  the  signal 
lies.  Here  "significant"  could,  for  example,  be  something  like  90%  or  95%  of  the  total 
energy.  The  precise  value  chosen  is  somewhat  arbitrary  and  will  usually  depend  on 
the  appHcation.  This  lack  of  a  unique  and  precise  definition  of  bandwidth  is  one  of  the 
reasorts  why  the  definition  of  broadband  is  also  imprecise.  By  using  this  definition  one 
can  write  the  fractional  bandwidth,  which  is  the  ratio  of  the  bandwidth  to  the  average 
frequency  within  that  band,  as 


A  = 


fu-fi 


m 


where  and  are  the  upper  and  lower  limits  of  the  specified  frequency  range 

respectively.  A  signal  is  considered  to  be  broadband  if  the  fractional  bandwidth  is 
greater  than  some  nominal  value,  e.g.,  A  >  0.1 ,  otherwise  it  is  narrowband.  A 
narrowband  signal  can  be  envisaged  as  one  that  approximates  a  sinusoid  as  A  0  . 

This  definition  of  a  broadband  signal  can  be  motivated  from  the  complex  signal  model 
which  assumes  that  the  signal  can  be  written  in  the  form 

^(0  =  Re{pW}  =  Re{a(0  exp(j2n/j)j  (2) 


where  a(t)  is  the  low  frequency  (complex)  envelope  and  the  exponential  term  contains 
the  high  frequency  carrier.  This  representation  is  valid  as  long  as  the  frequency 
spectrum  of  //(t)  does  not  have  any  negative  frequencies  (i.e.,  the  complex  signal  is  a 
good  approximation  to  the  so-called  analytic  signal,  see  appendix  A)  and  when  this  is 
die  case  the  signal  is  considered  to  be  narrowband.  Therefore  the  signal  is  broadband 
when  the  spectrum  of  ju(t)  contains  negative  frequencies,  i.e.,  Af  >  f.,  where 
A/  =  ^(/„  -  fi)  and  +  f,)  (see  Fig.  3).  This  means,  for  a  narrowband  signal, 

the  envelope  cannot  vary  at  a  rate  that  is  comparable  to  the  carrier  frequency  unless 
the  distortion  to  the  complex  signal  representation  is  deemed  acceptable.  In  other 
words,  for  broadband  signals,  the  carrier  frequency  term  cannot  be  separated  from  the 
envelope. 
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Figure  3.  Broadband  frequency  spectrum  of  a  complex  signal  with  negative  frequencies. 

2.2  Time-bandwidth  product  definition 

The  above  fractional  bandwidth  definition  of  a  broadband  signal  is  quite  general.  In 
active  sonar  a  more  specific  and  perhaps  useful  definition  can  be  formulated  in  terms 
of  the  signal  model  used  and  its  processing.  It  turns  out  tiiat  the  time-bandwidth 
product  of  the  received  signal  will  determine  whether  the  signal  should  be  treated  as 
broadband  or  narrowband.  To  see  this  consider  the  "slowly  fluctuating  point  target" 
model[2].  The  returned  signal  of  an  active  sonar  system  after  scattering  from  a  target 
can  be  modelled  as 


r{t)  =  Re{  a  -  T(i))  exp[727i/, (/  -  T{t))]  | ,  (3) 

where  a  is  a  complex  random  variable  describing  the  signal  attentuation  (caused  by 
transmission  loss  and  target  scattering)  and  the  range  imcertainty  of  the  target;  is 

the  carrier  frequency;  and  T(i)  is  the  time  dependent  roimd  trip  time  delay.  This  time 
delay  can  be  written  as  a  function  of  the  instantaneous  range  of  the  target,  R(t),  such 
that 


x{t)  =  ^F{t-\x{t)\,  (4) 

with  c  being  the  speed  of  sound  in  the  medium.  By  expanding  Eq.  (4)  in  a  Taylor 
series  about  the  initial  delay,  =  IR^jc,  and  assuming  that  the  target  moves  with 

constant  velocity,  then 

i-T(i)  =  i(i-To),  (5) 


where 
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C  +  V 

s  = - 


c-v 


(6) 


is  the  scaling  factor  and  v  =  v(to)  is  the  radial  velocity  of  the  target.  The  returned 
signal  of  Eq.  (3)  can  now  be  rewritten  as 

r{t)  =  Rej  a  exp[y27t(/,  +  fj){t  -  to)]  |, 

(7) 

where  the  Doppler  shift,  ,  was  introduced  through 


-fc=fc+fd  '^th 


fd  _  -2v/c 

fc  1  +  v/c  ■ 


(8) 


The  form  of  Eq.  (6)  shows  two  effects  on  the  returned  signal  which  are  caused  by  the 
target  velocity: 

•  The  signal  envelope  is  stretched  or  compressed  by  the  scaling  factor  s . 

•  The  carrier  frequency  is  Doppler  shifted  by  an  amount  . 

In  tihe  narrowband  approximation  the  scaling  of  die  signal  envelope  is  neglected.  To 
see  under  what  conditions  this  is  valid  compare  the  two  waveforms  a(i)  and  a{tls) . 
The  error  between  the  two  is  maximum  (see  Fig.  4)  when  t  =  T,  e.g.,  the  length  of  the 
pulse.  If  the  bandwidth  of  a(i)is  B  then  the  time  resolution  cell  size  is  1/5,  meaning 
that  the  signal  does  not  change  significantiy  during  this  time  interval.  Therefore  the 
error  between  a(T)  and  a(r/s)  is  small  if  [T  -  r/s|  « 1/5 .  This  can  be  rewritten  as 


TB« 


s 

5-1 


or 


TB«  — 
2v 


for  V  «  c . 


The  signal  is  considered  to  be  broadband  if  this  condition  is  violated. 


(9) 
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Figure  4.  (a)  the  difference  between  the  arguments  of  the  envelope  function  and 
(b)  the  resolution  cell  size  of  the  envelope  function. 


Note  that  the  condition  TB>\  may  not  be  a  sufficient  requirement  for  a  signal  to  be 
classed  as  broadband  in  accordance  with  the  above  definition.  This  can  be  illustrated 
by  the  following  example.  Consider  two  target  velocities:  (a)  50  knots  and  (b)  1  knot. 
By  using  Eq.  (9)  the  signal  wiU  be  narrowband  if  the  time-bandwidth  product  is  much 
less  than  30  for  case  (a)  and  750  for  (b).  If  the  required  time-bandwiddi  product  of  the 
signal  is  100  then  the  signal  is  broadband  for  case  (a)  whereas  it  can  be  considered 
narrowband  in  case  (b).  This  result  also  explains  why  radar  signals  with  substantial 
fractional  bandwidths  can  be  treated  in  the  narrowband  approximation  given  that  the 
propagation  speed  of  radar  signals  is  several  orders  of  magnitude  greater  than  that  of 
imderwater  soxmd. 


An  important  assumption  used  in  the  derivation  of  Eq.  (9)  is  that  the  target  maintains  a 
constant  velocity  during  the  observation  interval.  This  implies  that  narrowband 
processing,  which  is  limited  by  the  signal's  bandwidth,  can  be  used  to  model  a  system 
provided  its  accelerations  are  negligible.  Since  broadband  processing  is  required  when 
the  time-bandwidth  product  of  the  signal  is  of  the  order  of  c/2v,  the  performance  of  a 
narrowband  system  should  be  best  when  the  target  velocities  are  small.  Similarly, 
when  the  accelerations  of  the  system  become  non-negligible  the  performance  of 
broadband  processing  is  expected  to  fall.  This  limitation  of  broadband  processing  can 
be  expressed  quantitatively  as  follows.  The  acceleration  of  the  target  is  said  to  be  non- 
negligible  when  the  change  in  velocity  is  of  the  same  order  as  die  velocity,  i.e., 

Av  ~  V  w  —  —  .  Writing  the  acceleration  as  a  «  — ,  then 
2  TB  ^  T 


(10) 


It  is  evident  from  Eq.  (10)  that  large  accelerations  reduce  the  amount  of  processing 
time  available  for  broadband  systems. 
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In  summary,  velocities  degrade  the  performance  of  narrowband  systems  whereas 
accelerations  degrade  the  performance  of  broadband  systems.  The  advantage  of 
broadband  systems  is  that  they  allow  a  larger  processing  interval  which  results  in 
greater  gain,  better  noise  immunity  and  increased  range  resolution. 


3.  Broadband  Signal  Processing 

3.1  Matched  filtering 

The  matched  filter[3]  is  central  to  active  sonar  signal  processing.  Fundamentally,  the 
matched  filter  is  a  correlator  which  compares  the  received  signal  with  a  hypothesised 
signal  (or  set  of  signals).  The  output  of  the  matched  filter  gives  a  measure  of  how  well 
the  hypothesised  signal  (sometimes  referred  to  as  the  replica)  matches  the  received 
signal  as  a  fvmction  of  a  set  of  parameters,  usually  the  range  and  velocity  of  the  target. 
By  using  this  information  some  properties  of  the  target  can  be  inferred. 

More  specifically,  the  matched  filter  output  is  proportional  to 

x(t,  s)  =  y  * 


where 

x(t,s)  is  the  broadband  ambiguity  function[4], 

A:(i)  =  r(f)  +  «(if)is  the  received  signal  plus  noise  waveform, 

y(f)  is  the  hypothesised  signal  as  a  function  of  time  delay,  r ,  and  scale  factor,  s 
(see  Eqs.  (5)  and  (6)). 

Eq.  (11)  can  also  be  represented  diagrammatically  as  shown  in  Fig.  5. 


7 


DSTO-TR-0435 


x(t) 

received 


waveform 


>X(t,  s) 


Hypothesised 

replica 

waveform 


Figure  5.  Schematic  of  the  broadband  ambiguity  function 


If  the  received  signal  is  bandpass,  i.e.,  some  low  pass  function,  fip(f),  modulating  a 
carrier  with  frequency  f^. ,  then  the  received  signal  and  its  replica  can  be  written  as 

x{t)  =  Xjj.{t)QX^(j2'Kff]  and  y(t)  =  :>^i/>(^)exp(y27t//),  (12) 

so  that 

x(t,  s)  =  exp(^y27t/^  yip  *  exp(-  jlnfj)  dt  (13) 

The  magnitude  of  die  ambiguity  function  for  this  case  is  independent  of  the  carrier 
which  makes  it  easier  to  implement  digitally  since  only  low  pass  waveforms  have  to  be 
sampled.  Note  that  in  the  limit  as  5  — >  1 ,  the  ambiguity  function  in  Eq.  (13)  reduces  to 
the  well  known  narrowband  ambiguity  function  used  extensively  in  radar 
applications  [5]. 

The  ambiguity  function  is  a  useful  tool  for  analysis  and  provides  a  starting  point  for 
waveform  synthesis.  The  waveform  synthesis  problem,  which  consists  of  finding  the 
appropriate  signal  given  a  specified  ambiguity  function,  is  still  an  open  problem. 
Nevertheless,  from  the  experience  of  radar  signal  processing,  waveforms  can  be 
classified  into  three  broad  classes  of  ambiguity  functions:  the  ridge,  the  thumbtack  and 
the  bed  of  spikes.  The  ridge  function,  particularly  if  the  axis  of  the  ridge  is  oriented 
towards  the  Doppler  axis,  is  useful  for  the  detection  of  targets  with  unknown 
velocities.  The  technique  of  pulse  compression  is  one  means  of  realising  such  an 
ambiguity  function.  The  thmnbtack  function,  usually  considered  as  an  approximation 
to  the  "ideal"  ambiguity  function,  can  be  realised  through  the  use  of  pseudo-random 
waveforms.  Both  classes  of  ambiguity  functions  rely  on  waveforms  with  a  large  time- 
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bandwidth  product.  The  bed  of  spikes  function  is  obtained  whenever  a  train  of  pulses 
is  used.  Schematic  representations  of  these  ambiguity  functions  are  shown  in  Fig.  6. 


Figure  6.  Classes  of  ambiguity  functions:  (a)  the  ridge,  (b)  thumbtack  and  (c)  bed  of  spikes. 

3.2  Pulse  compression 

If  energy  considerations  are  ignored  then  the  short  pulse  has  a  number  of  desirable 
features  because  of  its  relatively  large  bandwidth.  These  include 

•  good  range  resolution;  the  ability  to  separate  multiple  targets  in  a  given  range 
interval  or  resolve  individual  scattering  centres  in  a  finite  sized  target 

•  good  range  accuracy,  which  follows  from  good  range  resolution 

•  mutipath  resolution 

•  a  shorter  minimum  range,  where  =  j  cT^  for  a  pulse  dxrration  of 

•  target  classification  through  better  feature  extraction 

•  Doppler  tolerance;  a  single  filter  matches  a  wide  range  of  Doppler  shifts 

The  use  of  explosive  charges  is  one  way  to  generate  short  pulses  with  adequate  energy 
though  this  may  not  be  always  operationally  desirable.  Alternatively  the  use  of  pulse 
compression  techniques  can  take  advantage  of  the  benefits  of  a  short  pulse  without  the 
constraints  of  peak  power  limitations. 

Pulse  compression[6]  occurs  if  a  waveform  with  a  nonlinear  phase  spectrum  is  passed 
through  a  filter  that  is  phased  matched  to  the  waveform.  A  signal  with  a  nonlinear 
phase  spectrum  is  one  whose  phase  is  a  nonlinear  function  of  time.  A  phased  matched 
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filter  is  one  whose  phase  response  is  equal  and  opposite  to  that  of  the  waveform.  (This 
is  just  a  special  case  of  the  matched  filter  described  in  the  previous  section.)  The  basic 
idea  is  best  illustrated  by  the  classic  example  of  linear  frequency  modulation  (LFM). 


VW\AP/ 

LFM  pulse  x(t) 


Figure  7.  Pulse  compression  using  a  tapped  delay  line  and  summation  configuration. 

Consider  the  setup  in  Fig.  7  where  the  tapped  delay  line  and  summation  configuration 
is  a  realization  of  a  matched  filter  for  the  input  LFM  waveform.  The  output  is  given  by 

y(r)  =  |;c(t)  h{t  -  x)  dx  (14) 

where  h(t')  =  x*  (-f)  is  the  matched  filter  condition  (for  simplicity  the  time  delay  for 
making  the  impulse  response  causal  has  been  neglected).  Let  the  input  LFM  waveform 
be  described  as 


^(0  =  expjy^TT 


for 


\,\<l 

‘  '  2 


(15) 


where  T  is  the  pulse  length  and  B  is  the  frequency  sweep  (this  is  a  good  measure  of 
the  signal's  bandwidth  if  >  50 ).  The  compressed  pvdse  is  then  given  by 


y{t)  =  Qxp{j2Kffi) 


r 

sin 

_ 1 

l  t)_ 

nBt 


for 


t\<T. 


(16) 


For  sufficiently  large  time-bandwidth  products  Eq.  (16)  has  zeros  for  f  «  —  for  \  t\<T , 

where  «  is  an  integer.  This  means  that  the  LFM  waveform  with  a  pulse  length  T  is 
compressed  to  a  pulse  whose  length  is  2/5 .  The  compression  ratio,  which  is  a  measmre 
of  the  amotmt  of  pulse  compression,  is  given  by  75/2,  proportional  to  the  time- 
bandwidth  product  of  the  uncompressed  signal. 
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From  the  above  analysis  it  is  evident  that  a  LFM  pulse  can  be  quite  long  in  duration. 
However,  since  each  instant  of  the  pulse  is  associated  with  a  unique  frequency,  the 
processor  can  combine  different  parts  of  the  pulse  in  such  a  way  so  that  the  output  is  a 
narrow  "spike"  which  accurately  specifies  the  target  delay  time. 

Another  form  of  pulse  compression  which  has  received  wide  application  is  the  phase- 
coded  waveform.  This  can  be  achieved  when  a  long  pulse  of  duration  T  is  divided 
into  N  subpulses  each  of  length  .  The  phase  of  each  subpulse  can  either  be  0  or  tc. 

If  these  two  possibilities  are  chosen  randomly  then  the  waveform  approximates  a 
noiselike  signal  with  a  thumbtack  type  ambiguity  function.  The  basic  idea  behind  the 
processing  is  that  the  sequence  of  short  pulses  is  combined  into  a  single  large  pulse 
with  many  of  the  small  pulses  as  sidelobes.  The  compressed  waveform  has  a  width 
and  an  amplitude  N  times  the  origmal.  A  special  phase-coded  sequence  is  the  Barker 
code  which  has  the  property  that  its  autocorrelation  fimction  has  equal  time- 
sidelobes[7]. 

It  should  be  noted  that  there  are  many  types  of  pulse  compression  waveforms.  These 
include  nonlinear  FM,  discrete  frequency  shift,  polyphase  codes,  compound  Barker 
codes  and  so  on.  AU  these  waveforms  have  their  particular  desirable  features  but  they 
come  at  the  cost  of  more  complicated  processing  and  the  practical  difficulties 
associated  with  their  implementation. 

3.3  Low  probability  of  interception 

A  necessary  condition  for  a  signal  to  have  a  low  probability  of  interception  (LPI)  is  that 
it  has  a  sufficiently  large  time-bandwidth  product.  Many  of  the  waveforms  used  for 
LPI  are  gimilar  to  the  pulse  compression  ones.  However,  to  be  effective  as  an  LPI 
signal,  the  design  of  the  signal  carrier  must  be  complex  enough  so  tiiat  its  timing, 
frequency  or  phase  pattern  cannot  be  predicted  by  an  unintended  interceptor.  Thus  the 
carrier  should  be  highly  variable  over  the  interval  in  which  the  message  or  baseband 
signal  is  nearly  constant.  This  means  that  the  modulated  signal  must  have  a 
bandwidth  larger  than  the  message  rate,  the  baseband  bandwidth.  In  general,  a  LPI 
signal  requires  its  total  bandwidtii,  5 ,  to  be  much  greater  than  the  message  rate  of 
1/r,  where  P  is  the  duration  of  the  baseband  signal,  i.e.,  TB»\.  The  objective  is  to 
make  the  broadband  signal  appear  noiselike  for  the  potential  interceptor.  The 
technique  used  to  achieve  LPI  signals  is  referred  to  as  spread  spectrum  modulation[8]. 

To  demonstrate  the  above  consider  a  signal  transmitted  using  some  sort  of  spread 
spectrum  modulation.  Assume  that  this  signal  at  the  intended  receiver  has  a  power,  5, 
with  bandwidth,  B ,  over  a  signalling  interval,  T .  Also  assume  that  tiie  interceptor  has 
a  bandwidth  IP  which  is  smaller  than  B .  Then  the  signal  to  noise  ratio  (SNR)  at  the 
interceptor  is 
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sm  = 

'  NjW  N,  BT 


(17) 


where  kj  is  a  factor  accounting  for  the  receiver  gain,  propagation  and  other  losses  at 
the  interceptor,  Nj  is  the  interceptor's  noise  spectral  density  and  the  intended  receiver 
power  can  be  written  as  S=  kEJT ,  with  E  being  the  transmitted  energy  over  the 
signaUing  interval  and  A  is  a  constant  of  proportionality.  Therefore  a  large  time- 
bandwidth  product  can  give  rise  to  a  relatively  small  interceptor  SNR. 


3.4  Beamforming 

The  discussion  so  far  has  assumed  that  the  received  signal  has  already  been  spatially 
filtered  and  is  in  a  form  suitable  for  temporal  processing.  The  objective  of  spatial 
filtering  or  beamforming  is  to  estimate  the  signal  coming  from  a  desired  direction.  It  is 
clear  from  the  preceding  discussions  that  broadband  signals,  whether  designed  for  the 
purpose  of  high  range  resolution  or  low  probability  of  interception,  are  characterised 
by  their  relatively  large  time-bandwidth  product  and  their  implications  and 
constraints  on  spatial  filtering  need  to  be  considered. 

Beamforming  is  based  on  the  idea  of  coherently  summing  the  received  signals  coming 
from  different  spatial  locations.  It  will  become  evident  from  the  following  discussion 
that  the  time  delay  between  received  signals  is  the  fundamental  quantity  in  broadband 
beamforming  rather  than  the  phase  shift  equivalent  which  is  used  extensively  for  the 
narrowband  case. 


sensors 


Figure  8.  Schematic  of  time  domain  beamforming. 
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To  illustrate  the  concepts  behind  the  beamforming  of  broadband  signals  consider  the 
linear  equi-spaced  array  depicted  in  Fig.  8.  For  simplicity,  assume  that  the  array  is  set 
up  to  coherently  sum  broadside  signals  (additional  pre-steering  delays  are  required  for 
the  coherent  summation  of  off-broadside  signals).  The  rectangular  pulses  shown  are  a 
generic  representation  of  the  envelope  of  the  broadband  signal.  The  actual  waveform 
could  be  anything  from  an  impulse-like  signal  to  long  duration  nonsinusoidal 
signals[9]  such  as  those  used  in  LPI  systems.  From  Fig.  8  note  that  all  the  sensors  are 
ensonified  as  long  as  0  =  sin"'  [cT^  jl] ,  where  is  the  pulse  duration  and  ^  is 

tihe  angle  of  arrival  of  the  wavefront  measured  from  the  broadside  direction.  As  6 
increases  beyond  6^  fewer  pxilses  will  get  summed  and  the  resultant  amplitude  drops 

until  01,  =  sin"' (cr^j/i)  where  T,^  is  the  pulse  repetition  interval.  At  this  stage  grating 
lobes  begin  to  make  their  appearance  as  pulses  from  the  next  period  get  summed. 

From  diese  simple  considerations  a  number  of  observations  can  be  made: 

•  The  angular  behaviour  of  the  beam  pattern  is  dependent  on  the  time  domain 
properties  of  the  waveform.  For  instance,  off  broadside,  the  resultant  waveform  at  the 
beamformer  output  becomes  increasingly  distorted  relative  to  the  broadside 
waveform. 

•  The  beamwidth  and  grating  lobes  depend  on  the  pulse  duration  and  the  pulse 
repetition  interval.  This  is  in  contrast  to  the  narrowband  formulation  of  beamforming 
where  the  beamwidth  and  grating  lobes  depend  on  the  carrier  frequency. 

•  The  shape  of  the  waveform  will  have  an  important  affect  on  the  performance  of  the 
beamformer.  It  was  found  in  Ref.  [10]  that  the  shape  of  the  signal  envelope  for  CW 
(continuous  wave)  pulses  is  an  important  factor  in  determining  the  directivity  pattern. 
Smooth  envelopes  give  a  better  directivity  than  ones  which  have  sharp  changes.  It  was 
also  fotmd  that  array  weighting  could  not  correct  for  badly  designed  signal  shapes  in 
this  respect. 

•  In  active  sonar  pulses  much  shorter  than  the  acoustic  aperture  are  often  used.  This 
can  result  in  the  whole  array  not  being  ensonified  simultaneously  and  hence  violating 
a  standard  assumption  used  for  narrowband  beamforming. 

•  The  use  of  time  limited  signals  in  active  sonar  suggest  that  it  is  more  pertinent  to 
describe  the  beam  pattern  in  terms  of  energy  distribution  with  respect  to  angle  rather 
than  power,  as  is  conventionally  done  in  the  narrowband  case.  Furthermore,  for 
pseudorandomly  coded  signals,  a  truer  performance  measure  of  the  receiver  response 
might  be  given  by  the  spatial  distribution  of  correlatable  energy[ll]. 

The  impression  given  by  the  foregoing  discussion  suggest  a  significant  departure  from 
narrowband  beamforming  concepts  and  techniques.  Nevertheless,  with  the 


13 


DSTO-TR-0435 


abovementioned  limitations  in  mind,  it  is  possible  to  generalize  the  single  frequency 
conventional  beamformer  to  one  that  can  process  time-limited  broadband  signals.  To 
see  this,  consider  a  linear  array  with  N  sensors.  The  beamformer  output  for  a  single 
frequency  component  (viz.,  the  narrowband  beamformer  case)  is 

2w.Z,.(/)exp(7(iTj)  (18) 

;=i 

where  Xi{/)  is  the  Fomier  transform  of  the  signal  at  the  zth  sensor,  is  the 
steering  time  delay  between  adjacent  sensors  and  w-  =  Wf  exp[ z  ,VU!)]  is  the 
sidelobe  shading  term  with  range  focussing  correction  (z,-  is  the  position  of  the  ith 
sensor,  X  is  the  wavelength  and  i?  is  the  range  of  the  target  measmed  from  some 
reference  sensor).  Then  the  beamformer  output  for  the  entire  frequency  spectnim  of 
the  signal  is 

N-\ 

12:  W,.  x{f)  exp(7(zT, ))  exp(y27r/i)  df .  (19) 

i=I 

This  formulation  of  a  broadband  beamformer  can  be  implemented  using  Fast  Fourier 
Transform  techniques[12]  as  shown  in  Fig.  9. 


Figure  9.  A  FFT  implementation  of  a  broadband  beamformer. 
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For  some  applications  requiring  the  transmission  and  reception  of  broadband  signals 
without  significant  deformation  of  its  frequency  spectrum  it  is  desirable  to  have  arrays 
capable  of  maintaining  a  constant  beamwidth  over  a  wide  frequency  range.  This  can 
be  achieved  by  configuring  the  array  such  that  each  sub-array  is  tuned  to  a  particular 
centre  frequency  so  that  the  desired  operating  frequency  range  is  covered  by  each  sub¬ 
array.  Examples  of  this  technique  can  be  fotmd  in  [13]. 


4.  Propagation  of  Broadband  Signals 


Up  to  now  the  discussion  of  broadband  signal  processing  has  implicitly  assumed  that 
broadband  signals  can  be  transmitted  through  the  medium,  i.e.,  the  underwater 
environment,  without  suffering  significant  distortion  to  their  frequency  spectra.  The 
purpose  of  this  section  is  to  see  under  what  conditions  this  assumption  is  valid.  More 
specifically,  to  determine  the  frequency  range  that  supports  broadband  propagation 
for  a  given  transmission  loss.  To  do  this  an  tmderstanding  of  how  signals  propagate 
through  the  medium  is  essential. 

From  a  study  of  ocean  acoustic  propagation  it  has  been  foimd  that  signal  degradation 
results  from  a  number  of  acoustic  loss  factors.  These  are  the  geometrical  spreading  loss 
(spherical,  cylindrical  etc),  the  volume  absorption  loss  (due  to  viscosity  and  chemical 
relaxation),  the  bottom  reflection  loss  (surface  reflection  loss  is  usually  considered  to  be 
negligible  due  the  large  impedance  mismatch  between  air  and  water)  and  the  stuface 
and  bottom  scattering  loss.  High  frequencies  are  foimd  to  be  strongly  affected  by 
boundary  scattering  and  the  microstructure  variability  of  die  ocean.  Therefore  volume 
and  scattering  losses  increase  with  increasing  frequency.  At  the  other  end  of  the 
spectrum  very  low  frequencies  are  found  to  be  strongly  affected  by  the  geoacoustic 
properties  of  the  sea  floor.  This  leaves  the  intermediate  frequencies  (of  die  order  of  10 
to  1000  Hz)  being  least  affected  by  losses  as  the  boundary  and  volume  effects  are 
minimized. 

In  order  to  get  a  quantitative  measure  of  die  "usable"  frequency  range  the  maximum 
and  minimum  detectable  frequencies  need  to  be  determined.  This  can  be  done  as 
follows: 

1.  Choose  a  sonar  system  configuration,  i.e.,  monostatic  or  bistatic,  the  type  of 
receiver  etc. 

2.  Determine  the  maximum  tolerable  transmission  loss  from  the  specifications  of  the 
sonar  system  via  the  sonar  equation[14]. 
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3.  Determine  the  transmission  loss  characteristics  of  the  medium,  e.g.,  the  ocean,  as  a 
hmction  of  frequency  and  range. 

By  combining  the  results  in  steps  2  and  3,  the  frequency  range(s)  with  tolerable 
propagation  losses  can  be  determined.  This  is  the  frequency  range  that  a  broadband 
signal  can  be  transmitted  without  suffering  significant  distortion  across  its  frequency 
spectrum. 

This  procedure  is  best  explained  through  an  illustrative  example.  Consider,  for 
simplicity,  a  monostatic  configuration  with  the  sonar  receiver  as  shown  in  Fig.  10. 


Figure  10.  Noncoherent  detector. 

The  signal  to  noise  ratio  of  the  receiver  determines  the  probability  of  detection 

/o(ax)  dx  (20) 


r 

1  jc  exp 

2 

for  a  given  probability  of  false  alarm,  Py^,  where  y  =  =  SNR(y) 

and/o(  )  is  the  zeroth  order  modified  Bessel  function.  Equation  (20)  can  be  presented 
as  a  set  of  receiver  operating  characteristic  curves  as  shown  in  Fig.  11. 
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Figure  11.  Noncoherent  detector  ROC  curves  for  false  alarm  probabilities  of  (a)  10-^  (b)  10^, 
(c)  10^  and  (d)  ICH. 


The  signal  to  noise  ratio  at  the  output  of  the  matched  filter  is  SNR{y)  =  lEJ Nq  ,  where 
£  is  the  received  signal  energy  and  jl  is  the  noise  spectral  density.  For  a  signal  of 
duration  T  and  bandwidth,  2B^ ,  at  the  beamformer  output,  the 


SNR{x)  = 


EIT  SNRjy) 
N,B,  ■  2BJ  ■ 


(21) 


By  specifying  the  required  and  then  SNR{y),  and  hence  SNR(x),  can  be 
determined.  To  continue,  asstune  the  following  set  of  system  parameters  in  a  noise 
limited  ocean  environment; 


1.  Pulse  length;  T=0.\s 
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2. 


Directivity  index:  D/  =  101ogjo 


f 

K 


10  dB,  where  the  size  of  the  effective 


aperture,  ,  of  the  receiving  array  is  taken  to  be  the  same  order  of  magnitude  as 
the  square  of  the  design  wavelength,  . 


3.  Source  level:  SL  =  200  dB  rel.  1  pPa  at  1  m. 

4.  Target  strength:  71S'=  10 

5.  Noise  level:  NL  =  50  dB  rel.  1  pPa  normalized  to  a  1  Hz  band. 

6.  Detection  threshold:  =  SNR(x)jg  =  SNR(y)j^ -10  logi^  T,  where  SNR(pc)^ 

has  been  normalized  to  a  1  Hz  band  to  be  consistent  with  the  noise  level 
specification. 


By  using  the  sonar  equation  the  maximum  tolerable  transmission  loss  can  be  written  as 


TL  =  UsL  +  TS-  NL„,  +  DI-DT^) . 


(22) 


Table  1,  by  using  Fig.  11,  shows  SNR{y)  for  a  given  set  of  and  The 

corresponding  set  of  transmission  losses  calculated  by  Eq.  (22)  using  the  above  set  of 
system  parameters  is  shown  in  table  2.  This  is  the  maximum  tolerable  propagation  loss 
for  a  given  specification  of  detection  and  false  alarm  probabilities. 


Table  1.  SNR  as  a  function  of  detection  and  false  alarm  probabilities 

SNR(y)dB  Pd  =  0.5  Pd  =  0.9 

Pfa  =  10^  12.375  14.5 

Pfa  =  10^  15.5  17.25 


Table  2.  Transmission  loss  as  a  function  of  detection  and  false  alarm  probabilities 

TLdB  Pd  =  0.5  Pd  =  0.9 

Pfa  =  1(H  73.8  72.8 

Pfa  =  10^  72.2  71.4 


1 

I 
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At  this  Stage  one  can  proceed  in  two  ways  to  get  the  range  of  frequencies  allowed  by 
the  limits  set  by  Eq.(22)  for  a  desired  minimum  range  of  detection.  The  first  is  to  use 
experimental  data  of  transmission  loss  which  is  a  function  of  frequency  and  range, 
such  as  that  shown  in  Fig.  12.  For  example,  one  such  experiment  was  performed  in 
coastal  waters  in  the  North  Atlantic  where  explosive  sources  were  used  to  map  the 
transmission  loss  as  a  function  of  range  and  frequency.  The  results  are  given  in  Fig. 
7(b)  of  Ref.  [15].  Using  these  results  for  a  maximum  tolerable  transmission  loss  of 
about  72  to  73  dB  gives  a  frequency  range  between  about  100  to  2000  Hz  for  a 
minimum  range  of  10  km.  This  is  for  a  source  and  receiver  depth  of  50  m  in  a  water 
column  approximately  100  m  deep. 


Figure  12.  Schematic  of  a  transmission  loss  contour  plot  as  a  function  of  range  and  frequency. 
The  actual  results  are  given  in  Ref  [15]. 


Alternatively,  if  experimental  data  is  unavailable,  the  transmission  loss  as  a  function  of 
range  and  frequency  can  be  obtained  by  modelling.  There  are  a  number  of  models 
available.  These  include  ray  theory  models,  range  independent  wave  models  (e.g., 
Fourier  Integral,  Normal  Modes)  and  range  dependent  wave  models  (e.g..  Coupled 
Modes,  Parabolic  Equation,  Finite  Element)  [16].  In  any  case,  aU  these  models  have  the 
common  objective  of  solving  the  wave  equation  with  a  harmonic  point  source  subject 
to  the  appropriate  boimdary  conditions: 


V''vj/(r,z)  + 


(0 


_c{r,z)} 


\\i(r,z)  =  - 


S(r)5(z-Zo) 


2nr 


(23) 


where  cylindrical  symmetry  has  been  assumed,  ip(r,z)  is  the  velocity  potential,  co  is 
the  angular  frequency  of  the  source  located  at  r  =  0  and  z  =  z^  and  c(r,z)  is  the  speed 
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of  soxmd.  For  the  purposes  of  continuing  tire  example  and  demonstrating  the  ideas 
involved,  a  simple  heuristic  model  will  be  used  instead  of  the  abovementioned  more 
sophisticated  models  which  are  computationally  intensive.  A  simple  transmission  loss 
model  is  given  by 

rZ,  =  10 log, o(r(,  xl0^)  +  101og,o(rxl0^)  +  fl(/)r  (24) 


where  is  the  transition  region  between  the  spherical  and  cylindrical  spreading  in 
km,  r  is  the  range  in  km  and  a{f)  is  the  attenuation  in  dB/km  accounting  for  the 
volume  and  scattering  losses  as  a  function  of  frequency.  To  get  an  estimate  of  the 
maximum  frequency,  take  H  (e.g.,  Zf  ~  100  m),  the  depdi  of  the  water  column  or 

duct,  and  a{f)  «  0.l(//ZZZz)^,  where  the  attenuation  is  assumed  to  be  dominated  by 
volume  effects  at  high  frequencies.  The  maximum  frequency  is  given  by 


?I-20-101og,.(r^xlO») 

■  ^kHz) 


(25) 


for  a  desired  minimum  range  .  For  a  minimum  range  of  10  km,  die  transmission 
loss  values  of  table  2  give  a  maximum  frequency  of  about  3.5  kFIz. 

To  get  the  minimum  frequency  note  that  sound  at  very  low  frequencies  tend  not  to  be 
trapped  in  a  duct  when  the  wavelength  becomes  comparable  to  the  dimensions  of  the 
duct  or  in  the  case  of  a  water  column,  interaction  with  the  bottom  becomes  significant. 
The  maximum  wavelength  is  given  by[17] 

yjn{z)-n{H)  dz  (26) 

where  n(z)  is  the  index  of  refraction  and  if  is  the  height  of  the  duct  or  water  column. 
For  a  constant  sound  velocity  profile 

or  (27) 

If  ZZ  —  lOOm  then  ~  100  Hz.  (It  should  be  noted  that  limitations  of  the  sonar 
hardware  could  set  a  higher  value  for  the  minimum  frequency.)  Therefore,  by  using 
the  simple  model,  the  useable  broadband  range  of  frequencies  is  expected  to  lie 
between  100  and  3500  Hz,  which  is  not  too  different  from  the  experimental  result. 

Although  the  above  example  only  used  a  minimum  range  of  10  km,  similar 
considerations  can  be  apphed  to  other  range  values  and  system  parameters. 
Qualitatively,  it  is  expected  that  the  detectable  range  of  frequencies  will  decrease  with 
increasing  minimtun  range.  Having  established  that  it  is  possible,  under  certain 
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conditions,  to  transmit  a  broadband  range  of  frequencies  without  serious  attenuation 
across  its  spectrum  the  next  step  would  be  to  determine  the  propagation  characteristics 
of  actual  broadband  waveforms. 

One  approach  is  to  treat  the  problem  of  broadband  propagation  as  the  same  as  that  of 
a  single  frequency  where  this  frequency  represents  the  geometric  mean  of  the  upper 
and  lower  limits  of  the  signal's  frequency  spectrum,  i.e.,  /„  =  V/i  fu  •  This  would  be  a 
reasonable  approximation  provided  the  signal's  bandwidth  is  not  too  large  and  the 
spectrum  remains  fairly  constant  over  the  bandwidth  of  interest. 

A  more  accurate  approach  would  be  to  use  existing  propagation  models  to  generate 
the  wave  solutions  for  each  frequency  of  interest  and  then  obtain  the  full  waveform 
solution  by  Fourier  synthesis[18],  i.e., 

f(r,  0  =  £  0  do)  (28) 

where  F(a))  is  the  temporal  spectrum  of  the  waveform  and  y/{r,z,(D)  is  the  spatial 
transfer  function  (in  terms  of  range  r,  depth  z  and  angular  frequency  (o)  obtained 
from  a  number  of  frequency  runs  of  a  wave  model  program.  Alternatively,  one  can 
attack  the  problem  directly  using  time  domain  methods[19].  It  shoxdd  be  noted  that 
both  these  methods  are  computationally  intensive. 


5.  Target  Classfication 


The  aim  of  target  classification  is  to  distinguish  one  target  from  another.  In  general, 
target  classification  requires  an  examination  of  the  details  of  the  echo  from  which  one 
can  hopefully  establish  a  signature  for  that  particular  target  in  question.  This  usually 
requires  a  larger  SNR  and  at  a  shorter  range  than  that  needed  for  detection.  There  are 
two  methods  of  target  classification  which  might  benefit  from  the  use  of  broadband 
signals.  They  are  resonance  scattering  and  high  resolution  range  profiling. 

The  resonance  scattering  method[20]  relies  on  the  target's  natural  frequencies  to  be 
excited.  A  complex  target  can  be  modelled  as  a  composite  of  standard  geometric 
shapes  (disks,  cylinders,  spheres  etc).  To  get  a  reliable  signature  of  the  target  it  is 
desirable  to  excite  as  many  resonances  as  possible.  Therefore  a  sufficiently  broadband 
signal  is  needed  to  estimate  most  of  the  target's  resonant  frequencies.  Once  a  set  of 
estimates  of  the  resonance  frequencies  has  been  obtained  then  a  template  can  be 
formed  for  future  classification  purposes.  Note  that  a  target' s  natural  resonances  are 
related  to  its  physical  size  so  the  largest  wavelength  that  can  be  expected  to  probe  the 
target  is  of  the  order  of  2  to  4  times  the  size  of  the  target.  On  the  other  hand  the  high 
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frequencies  will  excite  the  higher  modes  of  the  resonances  but  it  will  be  the 
fimdamental  resonant  frequency  which  will  give  the  greatest  response. 

The  range  profiling  method[21]  uses  high  range  resolution,  and  hence  large 
bandwidth,  signals  to  probe  the  target  so  that  all  the  major  scattering  centres  of  the 
target  can  be  resolved  individually.  This  creates  an  image  or  range  profile  of  the  target. 
From  this  range  profile  an  estimate  of  the  target's  size  can  be  made.  Further  target 
information  is  available  if  the  target  can  be  probed  at  multiple  aspect  angles. 

A  third  but  related  method  is  the  use  of  impulsive  sources  (explosions  or  short  pulses) 
to  determine  the  target  strength  over  a  wide  band  of  frequencies  [22].  This  has  the 
advantage  of  making  several  simultaneous  frequency  measurements  in  the  one 
operation.  By  treating  the  target  scattering  as  a  linear  process,  the  deconvolution  of  a 
traiismit-echo  pair  of  waveforms  gives  a  broadband  target  impulse  response  which 
characterizes  the  target.  Note  that  the  impulse  response,  or  its  transfer  function 
equivalent  in  the  frequency  domain,  is  also  a  function  of  target  aspect  angle  and  tilt 
angle.  So  a  more  complete  characterization  of  the  target  requires  impulse  response 
measxirements  to  be  made  over  a  number  of  these  angles.  Once  a  library  of  transfer 
functions  has  been  established  potential  targets  can  be  classified  accordingly.  In 
addition,  the  impulse  response  can  be  used  to  simiilate  target  echoes  for  any  given 
input  waveform  by  convolution. 


6.  Conclusions 


The  basic  motivation  for  using  broadband  sonar  is  to  improve  the  range  resolution 
capabilities  without  sacrificing  the  detection  performance,  which  is  governed  by  the 
amoimt  of  received  energy  from  the  target.  This  can  be  achieved  by  increasing  the 
bandwidth  of  the  signals  used.  From  this  improvement  in  range  resolution,  a  number 
of  other  desirable  features  foUow: 

•  better  range  accuracy 

•  improved  miiumum  range  of  detection 

•  greater  gains  in  the  signal  to  noise 

•  more  robust  estimates  of  target  parameters 

•  additional  information  for  target  classification 

•  possibility  of  multipath  resolution 

•  better  Doppler  tolerance 

•  lower  probability  of  signal  interception 

However,  all  these  features  come  at  a  cost  in  terms  of  increased  signal  processing 
complexity,  a  departure  from  the  narrowband  interpretation  of  measurements  and  the 
possible  limitations  in  the  maximum  detectable  range  and  useable  frequency  bands 
imposed  by  the  operating  environment. 
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Future  directions  for  broadband  active  sonar  research  are  expected  to  continue  in  the 
areas  of  signal  processing  where  improved  detection,  tracking  and  classification 
performance  in  highly  reverberant  environments  can  come  from  the  development  of 
tools  such  as  the  broadband  ambiguity  function,  wavelet  techniques  and  higher-order 
spectral  methods  (for  coloured  Gaussian  and  non-Gaussian  interferences). 
Implementation  of  beamforming  and  pulse  compression  techmques  wiU  continue  to 
attract  interest  with  the  emergence  of  extremely  fast  FFT  and  digital  signal  processing 
microtechnology.  Outside  of  the  signal  processing  domain,  work  is  needed  to  improve 
the  modelling  of  broadband  acoustic  propagation  in  terms  of  computational  speed  and 
accuracy.  In  closing,  it  is  clear  that  the  underlying  theme  which  connects  and  drives 
much  of  this  work  is  waveform  design. 
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8.  Appendix 


The  Analytic  Signal  Model 

Consider  the  following  signal  representation 

s{t)  =  Re{p(0}  =  RejflW  Qxp{j2nfj)} .  (Al) 

In  terms  of  its  Fourier  Transform,  S{f) ,  Eq.  (Al)  can  be  written  as 


s{t)=  [^S{f)exTp(j2nft)df 

=  j[  S{f)  Gxp{j2%ft)  df+^S{-  f)  exp(-  j2nft)  df  (A2) 

=  2  Rej  J[“  S{f )  exp(727r^)  d/| 

where  the  property,  S*{f)  =  S{-f)  for  a  real  signal  5(0,  was  used.  By  comparing 
Eqs.  (Al)  and  (A2),  the  complex  representation  is  given  by 

p(0  =  ^  2S[f )  exp{j2nft)  df  (A3) 

which  implies  that  the  spectrum  of  p(0  does  not  have  any  negative  frequencies.  If 
p(0  did  have  negative  frequencies  then  Re{p(0}  would  not  be  a  valid  representation 
of  s(t). 
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